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Abstract 

The method proposed by Pratt to derive recursion relations for 
systems of degenerate fermions [Phys. Rev. Lett. 84, 4255 (2000)] re- 
lies on diagrammatic techniques. This efficient formalism assumes no 
explicit two-body interactions, makes possible the inclusion of conser- 
vation laws and requires low computational time. In this brief report, 
we show that such recursion relations can be obtained from generating 
functions, without any restriction as concerns the number of conserva- 
tion laws (e.g. total energy or angular momentum). 

1 Introduction 

The statistical properties of canonical systems of fermions play a central 
role in the structural and radiative properties of multi-charged ions. Level 
densities, orbital populations, thermodynamic quantities, nuclear reaction 
rates and radiative opacities are crucial for the studies of laboratory and as- 
trophysical hot plasmas. Due to the large number of excited states at finite 
temperature, the calculation of the canonical partition function can become 
a difficult task, even for a reasonable number of particles. For instance, lean- 
ing on the fact that grand-canonical and canonical ensembles are equivalent 
in the thermodynamic limit, Kosov et al. [lj formulated recently a method 
to compute averages in the canonical ensemble based on calculations in the 
grand-canonical ensemble (which are much easier) . The technique presented 
by Pratt [2] to perform exact calculations of both canonical or microcanoni- 
cal quantities, involves the counting of all possible arrangements of fermions. 
It is based on recursive algorithms, which can be considered as extensions 
to techniques used for multiparticle symmetrization. Since the recursive 
method yields partition functions, it can be easily extended to the determi- 
nation of any statistical (average) quantity. The canonical partition function 
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Za-, for a system of A identical fermions populating a finite set of energy 
levels, ei, at an inverse temperature (3, reads: 

Z ^) = ^ E (si---SA\e^ H \V( Sl ---s A )), (1) 

' S 1 --:S A ,V{S) 

where the A distinguishable particles occupy the states si • • • sa, which 
are eigenstates of H, and where V(s) refers to the A\ permutations of 
si • • • s A - Pratt categorized the permutations diagrammatically, each dia- 
gram being expressed as a product of cycles. The sum over permutations 
is represented by the sum over diagrams. He extended the method to in- 
clude conserved quantities. For instance, conservation of the z component 
of angular momentum J z can be accomodated by summing over only those 
diagrams which yield a fixed sum: 

9 

M = J2 m u (2) 
i=l 

where rrii represent the momentum projection of state i (for instance, 
the number of states of a j orbital is g = 2j + 1). This problem of the 
determination of states resulting from the coupling of ./V particles can be 
carried out "directly" [H [3] , but the approach becomes tedious for systems 
with a large number of particles. The problem of listing the terms arising 
in a complex configuration can be solved from elementary group theory 
E] or by recursive techniques [7]. It was shown by Sunko and Svrtan 
[8j [9] and in a somewhat different way by Katriel [10\ lllj . that recursion 
relations can be obtained from a set of generating functions for the number 
of many-body sets, with the use of Gaussian polynomials (or q— binomial 
coefficients), In this brief report, we show that such a procedure can be 
extended at finite temperature without any restriction about the number 
of conservation laws. This enables one to recover the results of Pratt [2] 
and to draw a correspondence between the diagrammatic technique and the 
generating functions. The generating functions are introduced in Sec. [2j 
and the procedure to obtain the recursive relations is described in Sec. 
Section 0] is the conclusion. 

2 From q— binomial coefficients to generating func- 
tions 

It was shown by Sunko and Svrtan [5] that the multiplicities of the total 
angular momentum projections for n-particle systems are given by the co- 
efficients of some q— binomial coefficients [12] . Any polynomial in x may be 
written in terms of its roots X{. By collecting powers of x, one gets [T3] : 
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JJ(2!-2!j) = ^(-l) fc /fc(si,---,2; 9 j X 
1=1 k=0 

where fk is the k th elementary symmetric functions 



g-k 



(3) 



fk(xi, ■ ■ ■ ,x a ) = x h x---xx ik , (4) 

l<ii<— <ik<9 

where the summation extends over all ■ ■ ,ik} such that 1 < i\ < 
< i-k < 9- For instance, 



f3{xi,X 2 ,X 3 ,X4) = X±X 2 X 3 + X1X2X4 

+XIX3X4 + X 2 X 3 X4. 



(5) 



The number of states satisfying Eq. ([2]) is equal to the coefficient of 



(q t) in the polynomial [TTJ 



\{{l + q l t) = Y,t k fk(q,q\q\---,q 9 ) 

i=l 

which can be written 
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L 

r=0 



11(1 + 9**) = E8 



r(r+l)/2 



r=0 



(6) 



(7) 



where 



r 

L J q 



is the q— binomial coefficient (since we are focusing on 
fermions, g is an even integer): 

(1 - q9){l - q9- r ) ■ ■ ■ (1 - q9- r+1 ) 



J q 



(l-g)(l-g2)...(l- g r ) 



(8) 



The q— binomial coefficient is a generalization of the ordinary binomial 
coefficient, which corresponds to the particular case q = 1: 



0) 
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(0 


r 








q=l 





For instance, considering j = 7/2 and n = 3, we have 



-15/2 



g -15/2 + ,-13/2 + 2 g-U/2 + 3 ,-9/2 



+4 ? -V2 + 5g -5/J + 6g -3/2 + 6g -l/2 
+6 g 1 /2 +6 ,3/2 + 5g5 /2 +4g 7/2 

+3 ,f/* + 2q u ' 2 + q Vi / 2 + q 15 / 2 . 



(10) 
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3 From generating functions to recursive relations 



Following Pratt [2], let us consider a system of A identical fermions populat- 
ing a finite set of energy levels, ej, at an inverse temperature /?, and subject 
to the constraint M = J2i=i m i- l n that case, the partition function, noted 
Za,m{P)i is the coefficient of x M in the Taylor-series development of F A (x), 
defined as 



F A (x) 



1 d A F(x,z) 



A\ dz A 

where F(x, z) is the generating function 
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2=0' 



F(x,z) = ]ftl + 



i=i 



After a first derivative, one has: 



F A (x) 



1 



lvX' x ' 

qA-1 



dz A ~ 



2=0' 



which is equivalent to 



F A (x) 



x l e 
A-l rrf 



rrii -j3ei 



d A - 1 nLi(i + ^ 



-I3e kx m k ^ 



Q z A-\ (1 + ze -P^x m i) 



2 = 



(11) 



(12) 



(13) 



(14) 



Using Leibniz formula for the multiple derivative of a product of two 
functions, one obtains: 



F A {x) 



A\ 



E 



i=i 



A-l 

E 

71=0 
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A-l 



n 



d 



A-l-n 



Q z A-l-n (I + ze -i3e ix nir) 

Hil + ze-^x™*) 



k=l 



2=0 



(15) 



Since 
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one finds: 



d k 1 k\{-l) k a k 
dz 1 * (1 + az) ~ (l + az) k + 1 ' 



(16) 



F A - n {x), (17) 

which leads, since Za,m(P) is the coefficient of x M , to the following 
relation 

^A,M(/3) = ^EE e " n/3£l (- 1 ) n " 1 ^-^-^(/ 3 )' ( 18 ) 
n=l i=l 

which is exactly Eq. (6) of Ref. [2]. Therefore, the number of states 
with total energy E is obtained by (T — > oo, i.e. /3 — )• 0): 

A g 

Na,e,m = -7 E E^ -1 )™ -1 ^-™.^-™^^-™" 1 *- ( 19 ) 

n=l i=l 

Such a formalism can be extended to include additional constraints. For 
instance, in LS coupling, the total number of states with total energy E, 
reads : 

Na,e,m l ,m 8 = ^EEC-l)"- 1 * 

n=l i=l 

N A —n,E—nei,ML—nmn,Ms—nm s i ■ 

(20) 

Such recursion relations for the angular-momentum states are of great 
interest in detailed line-by-line calculations of atomic-structure and radiative 
opacity. For instance, they enable one to determine the size of the blocks of 
the hamiltonian matrix corresponding to a given value of the total angular 
momentum J, and to obtain the number of lines of a transition array. They 
can also play a role in the calculation of the rates of processes such as au- 
toionization. Applications related to atomic spectroscopy will be presented 
in a future paper. 

In the present work, as in Ref. [2j, no explicit two-body interactions 
are considered, i.e. it is assumed that the total energy of the system is a 
linear function of the numbers of particles. However, the interactions could 
be included, in an approximate way, either in a perturbative manner as sug- 
gested by Pratt [2], or following the variational Jensen- Feynman approach 
|14| . which consists in optimizing the one-particle energies in order to 
obtain the best agreement with the exact free energy. Such an approach 
is driven by the ability to obtain robust recursion relations |15|. [T6] for the 



F A {x) 



A 



]T(-ir-v 



-nfiei 



n=l 



x n.m,i 
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calculation of partition functions and introduces no further difficulties in the 
averaging process. 

4 Conclusion 

In this report, it was shown that the generating function method, formulated 
for instance through q— binomial coefficients, and which was introduced for 
the numbering of angular-momentum states of an arbitrary configuration, 
can be easily extended to calculate any average quantity in the canonical 
ensemble. The technique enables one to recover the recursive relations for 
the partition functions obtained by Pratt using diagrammatic techniques. 
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